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A brief review of modeling and simulation methods for a study of polymers at interfaces is pro- 
vided. When studying truly multiscale problems as provided by realistic polymer systems, coarse 
graining is practically unavoidable. In this process, degrees of freedom on smaller scales are elimi- 
nated to the favor of a model suitable for efficient study of the system behavior on larger length and 
time scales. We emphasize the need to distinguish between dynamic and static properties regarding 
the model validation. A model which accurately reproduces static properties may fail completely 
when it comes to the dynamic behavior of the system. Furthermore, we comment on the use of 
Monte Carlo method in polymer science as compared to molecular dynamics simulations. Using the 
latter approach, we also discuss results of recent computer simulations on the properties of polymers 
close to solid substrates. This includes both generic features (as also observed in the case of simpler 
molecular models) as well as polymer specific properties. The predictive power of computer simu- 
lations is highlighted by providing experimental evidence for these observations. Some important 
implications of these results for an understanding of mechanical properties of thin polymer films 
and coatings are also worked out. 



I. INTRODUCTION 

Due to their ubiquitous presence in everyday life, the 
importance of polymers can hardly be over estimated. 
Indeed, polymer science has had a major impact on the 
way we live. Just 50 years ago, materials we now take 
for granted did not exist at all. Due to their structural 
complexity, polymers are generally not crystalline at low 
temperatures. Rather, they exhibit an amorphous, glassy 
structure. The glass transition (the transition from a liq- 
uid to an amorphous solid) is thus an important concept 
when it comes to an understanding of the properties of 
polymer systems. 

As to the application, polymers are often used as pro- 
tective coatings in many fields ranging from the car in- 
dustry (corrosion resistance) to microelectronics (protec- 
tion against thermal as well as mechanical load) [H, 0, 0] . 

In such situations, the polymer is always in contact 
with a solid substrate or even confined between two solid 
plates (film geometry). An important piece of informa- 
tion for materials design is therefore how the thermal and 
mechanical properties of a polymer system are affected 
by interaction with a solid substrate. 

In addition to its technological importance, the investi- 
gation of polymers close to solid surfaces is also of great 
theoretical interest. One is, for example interested to 
understand how the broken symmetry of the space in 
the proximity of the substrate as well as related ener- 
getic issues modify static and dynamic properties (chain 
conformation, diffusion, reorientation dynamics, etc.) of 
polymers close to the interface. A particularly interest- 
ing aspect is also how the glass transition is affected by 
the presence of the solid substrates. 

However, the presence of a wide range of time and 
length scales makes the study ofpolymer systems a very 
challenging multiscale problem jj]. In the simplest case 



of linear homopolymers, each macromolecule contains N p 
identical repeat units (monomers), connected to form a 
chain. In experiments, the chain length may vary in the 
range 10 2 < N p < 10 6 . This implies that the average 
size of a polymer, measured for instance by the radius 
of gyration R g [!, 0], varies between R g ~ lOd up to 
R g ~ lOOOd, where d denotes the diameter of a single 
monomer (d being of the order of a few A). 

These different length scales are reflected in the par- 
ticular features of a polymer melt. In the melt the 
monomers pack densely, leading to an amorphous short- 
range order on a local scale and to an overall low com- 
pressibility of the melt. Both features are characteristic 
of the liquid state. Qualitatively, the collective struc- 
ture of the melt thus agrees with that of non-polymeric 
liquids. Additional features, however, occur if one con- 
siders the scale of a chain. A long polymer in a (three- 
dimensional) melt is not a compact, but a self-similar 
object @, 0, Q. 

It possesses a fractal 'open' structure 
which allows other chains to penetrate into the volume 
defined by its radius of gyration. On average, a polymer 
interacts with y/N p other chains, a huge number in the 
large- N p limit. This strong intcrpenctration of the chains 
has important consequences. For instance, intrachain ex- 
cluded volume interactions, which would swell the poly- 
mer in dilute solution, are screened by neighboring chains 
IS 0, 0, 0, EH E3 > although nontrivial nongaussian ori- 
entational correlations due to excluded volume remain 




A polymer in a melt thus in many respects behaves 
on large scales as if it were a random coil, implying that 
its radius of gyration scales with chain length like R g ~ 
\/N p . Furthermore, the interpenetration of the chains 
creates a temporary network of topological constraints 
@, 0, 0, [HI- These entanglements greatly slow down 
the chain dynamics and render the melt in general very 
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viscous compared to low-molecular weight liquids. 

In this article, we will provide a brief overview on com- 
puter simulations of polymers close to interfaces, while 
at the same time briefly touching upon some of the 
techniques aimed at a reduction in the computational 
cost related to the presence of multiple length and time 
scales. It should, however, be emphasized here that the 
present survey is far from being exhaustive, neither in 
regard to the modeling of polymer melts nor in regard to 
specific computational aspects. We refer the interested 
reader to extended reviews on these topics, e.g. refer- 
ences p EI EE EE El El El ■ 

In principle, a full multiscale simulation approach con- 
structs the connection between the quantum description 
of interactions (Fig. 1) over the all- atom description of a 
chemically realistic model but using effective classical po- 
tentials up to the coarse-grained level of mesoscale mod- 
els explicitly. Such an explicit multiscale approach which 
derives effective interactions on the mesoscale ab initio 
still is an ambitious and difficult task [IE HH, HE H3 and 
most of the work in the literature considers only partial 
steps of the full problem. Thus, in this article, we will 
concentrate on the description in terms of coarse-grained 
bead-spring type models lacking explicit chemical details, 
and we shall describe the steps (Fig. 1) to derive these 
models in a schematic way only. 



II. COARSE-GRAINING 

On a fundamental level, interaction forces originate 
from the adaption of the electronic degrees of freedom 
to the positions of the nuclei. It may therefore appear 
natural to model polymer melts via the Car-Parrinello 
method [24|. This method is a molecular dynamics (MD) 
technique |25l. l2u| which allows the electrons to adiabati- 
cally follow the motion of the nuclei, thereby replicating 
the energy landscape that the nuclei feel at any instant 
of their motion. Using recent extensions of the method, 
one is now able to study system sizes up to about 1000 
nuclei for about 10 ns [27|. This time, however, barely 
suffices to equilibrate the system at high temperature in 
the liquid state [IE EE HI- 

Some kind of coarse-graining procedure is therefore 
necessary in order to adequately describe statistical me- 
chanical properties of polymer systems on time and 
length scales of interest. Such a procedure usually con- 
sists of elimination of fast degrees of freedom by incorpo- 
rating them in effective potentials [3, HE Ell ■ 



A. Atomistic models 

A first degree of simplification may consist in replac- 
ing the electronic degrees of freedom by empirical poten- 
tials for the bond lengths, the bond angles, the torsional 
angles and the nonbonded interactions between distant 
monomers along the chain ('quantum level — > atomistic 
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FIG. 1: A schematic view of different scales which can be 
focused upon when describing polymers. The quantum level 
takes account of the electrons to calculate the interactions 
between the nuclei. On the computationally less demanding 
atomistic level, the electronic degrees of freedom are replaced 
by a force field [29[. Computationally still less demanding 
than atomistic models are simulations at the coarse-grained 
level. Here, a monomer is associated with a spherical site and 
the realistic potentials are replaced by simpler ones [ID, H3 . 
0,1. 



level', see Fig.[T]). This step introduces a 'force field', i.e., 
the form of the potentials is postulated and the corre- 
sponding parameters (e.g. equilibrium bond length, force 
constants, etc.) are determined from quantum-chemical 
calculations and experiments [IE HE] ■ 

Several such force fields have been proposed through- 
out the past decades for both explicit atom models and 
united atom models. An explicit atom model treats ev- 
ery atom as a separate interaction site, whereas a united 
atom model l ump s a small number of real atoms together 
into one site [3 HE H3] ■ Typical united atoms are CH, 
CH2, and CH3. The reduction offeree centers translates 
into the computational advantage of allowing longer sim- 
ulation times. With a time step of ~10~ 15 s — compared 
to ^10 -17 s for the Car-Parrinello method — a few thou- 
sand united atoms can be simulated over a time lapse of 
several /is, about an order of magnitude longer than an 
explicit atom simulation of comparable system size. 

Both explicit atom models and united atom models 
have been used in the study of glass-forming polymers 
(see e.g. [HI, E3 for reviews on older work). Current 
examples include polyisoprene (explicit atom; [33|, E3), 
atactic polystyrene (united atom; [IE HE H3) an d c * s_ 
tran s 1 . 4-pol ybu tadiene (united and explicit atom mod- 
els; [IE SlHE 5E EE El)- Certainly, the ultimate 
objective of these modeling efforts is that the simulation 
results lend themselves to a quantitative comparison with 
experiments. Such a comparison may, however, require 
a careful fine-tuning of the force field. For the family 
of neutral hydrocarbon polymers the optimization of the 
torsional potential appears particularly crucial. Not only 
the position and the relative depth of the minima, but 
also the barriers between them should be accurately de- 
termined, as local relaxation processes, involving tran- 
sitions between the minima, are exponentially sensitive 
to them. In extreme cases, imprecise barrier heights may 
seriously affect the dynamics while leaving structural fea- 
tures of the melt unaltered. 

Such an example is shown in Fig. [5] The figure com- 
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FIG. 2: Simulation results for cis-trans 1,4-polybutadiene 
(adapted from reference [3!|). Main panel: Collective static 
structure factor, S(q), and single-chain structure factor, w(q), 
versus the modulus of the wave vector q at T = 273 K. 
Two united atom models are compared: a chemically real- 
istic model (CRC) and the same model but without torsional 
potential (FRC). The vertical arrows indicate the q- values 
associated with the radius of gyration and with the first max- 
imum of S(q) ('amorphous halo'). The maximum occurs at 
q* ~ 1.47 A -1 . In real space, this value would correspond 
to an intermonomer distance of « 4.3 A which is roughly 
compatible with the average Lennard-Jones diameter of the 
model (d ~ 3.8 A). Inset: Mean-square displacement go(t), 
averaged over all monomers, versus time for the CRC and 
FRC models at T = 273 K. The horizontal line indicates the 
radius of gyration R g 2 = 218 A 2 (which is found to be the 
same for both models [40lp. 



pares simulation results for two models of a polybutadi- 
ene melt: a carefully validated united atom model which 
reproduces the experimentally found structure and dy- 
namics of the melt, and the same model with the tor- 
sional potential switched off. Apparently, suppression of 
the torsional potential has no influence on the structure, 
but considerably accelerates the monomer dynamics. 

The above example demonstrates that different poten- 
tials may lead to realistic representations of structural 
properties, but to diverging predictions for the dynam- 
ics. Such an observation is not limited to polymers, but 
was also made e.g. for amorphous SiO-2 [43]. It suggests 
that the design of a chemically realistic model, aiming 
at a parameter-free comparison between simulation and 
experiment, should involve information about both the 
structural and the dynamic properties. 



B. Generic models 

In many cases, one is interested in a basic understand- 
ing of polymeric features rather than detailed description 
of a specific model. For uncharged linear polymers, these 
features are presumed to be the dependence of mate- 



rial properties on chain connectivity, excluded- volume in- 
teractions, monomer-monomer attractions and/or some 
stiffness along the chain backbone. In this case, it ap- 
pears permissible to forgo fast degrees of freedom (bond 
length and bond angle vibrations, etc.) in favor of a 
coarse-grained model. A special type of coarse-grained 
models are so-called 'generic models' [l4[- Various such 
generic models have been studied in the literature (for 
reviews see references [13, 0, EE EH ) ■ Due to the lack 
of space, in the following we present only one of these 
models in more detail, which was used in recent Simula- 
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53, 54, 55, 56, 57, 58, 59, 6 



63L 1 641 . |65| . Certainly, this choice is biased by our 
own experience. 

Indeed, the issue of polymers at interfaces in par- 
ticular and confined liquids in general has received 
considerable attention in simulation studies. Various 
systems — simple liquids jr36l Ir37l l68l.[69l . [70l . l7lll . hydrogen- 
bonded or molecular liquids [72, u3 , silica 



[75j,[76j,[77|,[78j,[79j,[8i]|,[81 



t ■? ■ '? ■ ■? ■ '? ■ ■? ■ — ■> | i -nam- 
ing geometries — pores [66J, [72 




74j . polymers 
— and confin- 
[73j jfillers in g lass-for ming 
thin films J67|, M, [Zfl IzE M, M, 
3, 0, ^9, ^91 — have been consid- 



matrice s I69L 176 . 177 

ered. References provided here as well as throughout 
this paper are hoped to, at least partially, compensate 
this shortcoming. 



III. A BEAD-SPRING MODEL FOR POLYMER 
MELTS 

In 1990 Kremer and Grest proposed a versatile 
bead-spring model for the simulation of polymer systems. 
The 'Kremer-Grest' model has ever since been deployed 
to investigate numerous problems in polymer physics, in- 
cludin g re laxation processes in polymer solutions [88[ and 
mel ts |89l l9(il [9l| or the behavior of polymer brushes 
[92l |93|. to name just a few. 

In a variant of this model, proposed by Bennemann and 
coworkers [47| , the chains contain N p identical monomers 
of mass m. All monomers, bonded and nonbonded ones, 
interact by a truncated Lennard-Jones (LJ) potential 



4e [{d/r) 12 




{d/rf 



for r < 7' c 
else . 



(1) 



The parameter r c = 2r m j n is the cut-off distance, where 
i"mm = 2 1 ' 6 d is the minimum of Eq. ([T]). In contrast to 
the original version of the Kremer Grest model, where 
r c = ?*min (leading to purely repulsive intcrmolccular in- 
teractions), the cutoff distance proposed by Bennemann 
is motivated by the wish to work with a potential that is 
as short-ranged as possible while yet including the major 
part of the attractive van-der-Waals interaction. Even 
though attractive interactions are not expected to appre- 
ciably affect the local structure in a dense melt, they have 
a significant effect on thermodynamic properties. Fur- 
thermore, they arc important for simulations of e.g. the 
phase behavior of polymer solutions 0, [95| , thin films 
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with a film-air interface [96l . |97| or crazing in polymer 
glasses HHH. 

The chain connectivity is ensured by a FENE (finitely 
extensible non-linear elastic) potential 



£/fene(0 = -ifci?ohi 



1-f- 



Rq = 1.5d, k = ^ 
d z 

(2) 

The FENE potential diverges logarithmically in the limit 
of r — > Rq ('finite extensibility') and vanishes paraboli- 
cally as r — > ('elastic behavior'). The superposition of 
the FENE- and the LJ-potcntials yields a steep effective 
bond potential with a minimum at r\> w 0.96c? (see for 
instance Ref. [oOl]). The difference between rb and r m ; n 
is crucial for the ability of the model to bypass crystal- 
lization and exhibit glass-like freezing in. 



A. Approximate mapping to real units 

The parameters of Eq. ([I]) define the characteristic 
scales of the melt: e the energy scale, d the length scale, 
and tlj = (mcP/e) 1 / 2 the time scale. In the following, 
we utilize LJ-units. That is, e = 1, d = 1, and m = 1. 
Furthermore, temperature is measured in units of e/fce 
with the Boltzmann constant fee = 1- 

Although reduced units are commonly e mployed in 
simulations and are of technical advantage [lOCt Il0l| . 
it might still be interesting to obtain a feeling how they 
translate into physical units. Such a mapping to real sys- 
tems has recently been carried out by Virnau [§5|, Il02l | 
and by Paul and Smith (29| . Virnau explored the phase 
separation kinetics of a mixture of hexadecane (C16H34) 
and carbon dioxide (CO2). By identifying the critical 
point of the liquid-gas transition in hexadecane with that 
of bead-spring chains containing 5 monomers they found 
d ~ 4.5 x 1CT 10 m and e ~ 5.8 x KT 21 J. Paul and 
Smith used the data on the dynamics of chemically real- 
istic models for noncntangled melts of polyethylene and 
polybutadiene and obtained tlj ~ 0.21 x 10~ 12 s. These 
values for d, e, and tlj are compatible with the estimates 
obtained by Krcmcr and Grcst when comparing the dy- 
namics of entan gled bead-spring melts to real polymers 
(see table III of [87|). 



B. Choice of the chain length 

In polymer glass simulations the chain length is usually 
chosen as a compromise between two opposing wishes: 
On the one hand, it should be sufficiently large to sep- 
arate the scales of the monomer and of the chain size 
so that polymer-specific effects (or at least the onset 
thereof) become observable. On the other hand, compu- 
tational expedience suggests to work with short chains. 
Because the simulations aim at following the increase of 
the monomeric relaxation time To with decreasing tem- 
perature over as many decades as possible, slow relax- 
ation processes, already present at high temperatures (T) 




FIG. 3: Snapshot of a polymer system between two substrates 
of triangular lattice structure (only 40 chains out of 200, each 
containing N p = 10 monomers, are shown). 



due to entanglements, should be avoided. Thus, the chain 
length should be smaller (or at least not much larger) 
than the entanglement length N c . Extensive studies of 
the Kremer-Grest model show that N c « 32 Shorter 
chains exhibit Rouse-like dynamics 



r(N p ) = r A~ 



(3) 



As the Bennemann model is expected to have a similar 
N c , the chain length N p = 10 was proposed as a possi- 



ble compromise [47] 
subse quen t studies 1 
melts MiSi 



This chain length was used in all 
iertainin g to glass-forming polymer 

I, HH, MM M HE HI, M M M 



C. Including solid substrates 

Even though real substrates can have a complex struc- 
ture, it appears natural in the spirit of the polymer mod- 
els discussed above also to treat the substrate at a generic 
level. One obvious feature is its impenetrability. So a 
minimal model must at least respect monomer-substrate 
excluded volume interactions. Further generic features 
could be some surface roughness and adhesive power. 
Based on this reasoni ng, s imul ations often model the 
substrate as a crystal |l03l Il04| made of particles that 
interact with each other and with the monomers by LJ- 
potcntials. 

Such crystalline substrates may be implemented by 
tethering the substrate ('wall') atoms to the sites of a 
triangu lar lattice [oil ] via harmonic springs ('Tomlinson 
model' pMHo^l 



Ut(v) = h T (r- r cq ) 2 , k T = 100 (LJ units). (4) 



Here, r oq denotes the equilibrium position of an atom 
on the triangular lattice and fcx the spring constant. The 
substrate atoms are LJ-particlcs that interact with each 
other and with the monomers. The parameters (e and d) 
for these interactions — wall-wall and monomer-wall — are 
the same as in Eq. (JTJ) . 

Since the spatial arrangement of the substrate atoms 
may have a strong influence on the properties of the melt 
in the very vicinity of the substrate, it is interesting to 
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also study the effect of substrates with an amorphous 
structure. This can be achieved in a way quite similar to 
the implementation of crystalline walls (Fig.[3|); the only 
difference being the random (instead of regular) distri- 
bution of substrate atoms. 

If one is only interested in the average force which the 
substrate exerts on a monomer, one may treat the wall 
as a continuum and integrate over the parallel (x, in- 
directions and the vertical one up to the wall-melt inter- 
face. Carrying out this calculation for the LJ-potential 
one obtains 



U w (z) 



d\ 3 ' 



(5) 



where e w denotes the monomer-wall interaction energy 
and / w is a constant. While the second attractive term 
is im portant if one wants to study po lymer adsorption 
1 1 07t ] or wetting phenomena [l08l Il09| | , the first term of 
Eq. §5$) suffices to impose a geometric confinement. This 
is the stance we have adopted in most of the simulations 
on supercooled polymer films [13, [H, [H, [62J . 



IV. MOLECULAR DYNAMICS VERSUS 
MONTE CARLO 

In the framework of computer simulations, it appears 
natural to address dynamical problems via MD tech- 
niques. However, if one is interested in equilibrating 
long-chain glass-forming polymer melts at low T, MD 
might not be the most efficient approach. The realistic 
molecular dynamics has the drawback that the available 
simulation time is often not sufficiently long for an equili- 
bration of the polymer configuration at low temperatures 
or for long chain lengths. 

Therefore, one migh t envisage resortin g to Monte 
Carlo (MC) techniques [HE IlHlIll E3- The strate- 
gic advantage offered by this method is the number of 
ways in which MC moves may be designed to explore 
configuration space. The hope is to find an algorithm 
that, freed of the need to capture the real dynamics, ef- 
ficiently decorrelates the configurations of glass-forming 
polymer melts at low T. This demand on the algorithm 
appears to exclude the simplest MC technique, the ap- 
plication of only local MC moves, as a possible candi- 
date. A local MC move consists of selecting a monomer 
at random and in attempting to displace it by a small 
amount in a randomly chosen direction [45|. Not only 
should the local character of coordinate updating share 
the essential problematic features of the (local) molecular 
dynamics at low T or for large N p , but also may one ex- 
pect that local MC moves will yield an unfavorably large 
prefactor of the relaxation time due to their stochastic 
character. This conjecture is base d on an observation 
made by Gleim and coworkers |l!4| . They compared the 
relaxation dynamics of a glass-forming binary mixture 
simulated, on the one hand, by MD and, on the other 
hand, by a stochastic (Brownian) dynamics (which is in 



some respect similar to MC). They demonstrated that, 
although the structural relaxation at long times is the 
same for both methods, MD is roughly an order of mag- 
nitude faster than the stochastic dynamics. 

However, MC moves need not be local. They can be 
tailored to alter large portions of a chain. A prominent 
example of such nonlocal moves is the confi guration-bias 
Monte Carlo (CBMC) technique [H,[lOl|. Application 
of this technique to dense polymer systems in the canon- 
ical ensemble usually involves the attempt to remove a 
portion of a chain starting from one of its monomers that 
is randomly chosen and to rcgrow the removed portion 
subject to the constraints imposed by the local potential 
energy If successful, this implies a large modification of 
the chain configuration, thereby promising efficient equi- 
libration. However, Bennemann found that even in the 
limit where only the end is reconstructed ('smart repta- 
tion'), CBMC is inferior to ordinary MD [43]. In a dense 
melt, the probability of inserting a monomer becomes 
vanishingly small anywhere except at the position where 
it was removed. So, the old configuration of the chain 
is just restored. This trapping of the chain makes the 
relaxation become very slow. 

Thus, successful nonlocal chain updates in dense sys- 
tems should involve moves that do not require (much) 
empty space. A promising candidate is double-bridging 
algorithms which were succe ssfully em ployed in simula- 
tions of polye thyle ne chains jl!5l Ill6| . of the Kremer- 
Grest model 117L and of a lattice model, the bond- 
fluctuation model The basic idea of the algorithm 
is to find pairs of neighboring chains which one can de- 
compose into two halves and reconnect in a way that 
preserves the monodispcrsity of the polymers. Such a 
connectivity-altering move drastically modifies the con- 
formation of the two chains and thus strongly reduces 
the slowing of the dynamics due to large values of N p . 
However, if we attempt to repeat this move over and over 
again on the melt configuration we started with, a suc- 
cessful double-bridging event is likely to annihilate one of 
its predecessors by performing the transition between two 
chains in the reverse direction. To avoid this inefficiency 
this nonlocal chain updating should be complemented by 
a move which efficiently mixes up the local structure of 
the melt. At low T, efficient relaxation of the liquid struc- 
ture calls for a method which alleviates the glassy slowing 
down in general. Thus, any algorithm achieving this aim 
in non-polymeric liquids should also accelerate the equi- 
libration of glassy polymer melts, provided that it can be 
generalized to respect chain connectivity. At present, no 
technique has been established to fully solve this problem 
(see Ref. jl!8l | for a topical review). However, promising 
candidates a ppear to be 'parallel tempering' [1 19l ll2C| 
(see however |l2l| ). a recent MC approach proposed by 
Harmandar is and Th eodorou combining scission and fu- 
sion m oves Il22l Il23l ] or variants of 'Wang-Landau sam- 
pling' pTll . ll24| . 
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FIG. 4: Molecular dynamics simulations of substrate effects 
on the structure (a) and dynamics (b) of a model polymer 
melt described by Eqs. ([T]) and (J2)). (a): Density profile (nor- 
malized to the liquid density at infinite distance) versus dis- 
tance z from the substrate, (b): Local mean square displace- 
ments (MSD) versus time at two different distances, z, from 
a solid substrate Ugl. 
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FIG. 5: Growth of both the strength and the range of sub- 
strate effects upon cooling, (a): Relaxation time versus dis- 
tance from the substrate for a model polymer melt described 
by Eqs. ([T]) and ([2]) at various temperatures (Lennard-Jones 
units). The horizontal dashed lines indicate the bulk val- 
ues (expected at large distances from the substrate), (b): 
Monomer number density profile for the same range of tem- 
peratures [5g| | . 



V. POLYMERS AT INTERFACES: SOME 
SALIENT FEATURES 

In order to demonstrate the capabilities of the simple 
model described above, we focus here on our recent sim- 
ulation studies. Compared to the earlier Monte Carlo 
studies of the bond- fluct uatio n lattice model [Ifl, re- 
viewed in references [3, Il25l |. our work deployed MD 
simulations to explore the features of a continuum model, 
spatially confined to a slab geometry [5^, [|5(| [Ml, HI] ■ 

As an example, Fig. |4] compiles results of molecular 
dynamics simulations of the polymer model of Eqs. ([1]) 
and @ in the proximity of a solid substrate (at distances 
of a few monomer diameters), showing the strong het- 
erogeneity induced by the substrate both in the static 
structure (exemplified by the local density) as well as in 
transport properties of the liquid. The figure also illus- 
trates the influence of different solid structures such as 
a perfectly smooth substrate as compared to substrates 
with atomic scale corrugation (amorphous as well as crys- 
talline) . 



As shown in Fig. 01 the liquid density exhibits oscilla- 
tions in the proximity of a substrate. These oscillations 
are of comparable magnitude both for an ideally flat and 
an amorphous corrugated substrate. Despite this simi- 
larity in the behavior of the local density, the effect on 
dynamics is of opposite nature. While a perfectly smooth 
substrate (no corrugation at all) leads to an acceleration 
of diffusive dynamics, the dynamics of monomers close 
to an amorphous substrate is slowed down. A qualitative 
understanding of this, at first glance unexpected, behav- 
ior can be obtained by invoking the concept of effective 
friction gHH. 

Figure [4] also underlines the fact that the spatial ar- 
rangement of substrate atoms may play a crucial role in 
the properties of the adjacent liquid. Strong layering of 
liquid particles has a dynamic counterpart manifest in 
a temporary arrest within the liquid layer closest to the 
substrate (note the extended plateau in the mean square 
displacement). 

Substrate effects both on the packing behavior as well 
as on dynamics can become quite dramatic when cool- 
ing the liquid towards the freezing temperature. This 
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FIG. 6: Scaling plot of T c (h) and T e (h). Here h is the film 
thickness and T c (h) is the so called critical temperature of 
the mode coupling theory determined in MD simulations of 
the polymer model specified via Eqs. {1} and ((2| confined be- 
tween two purely repulsive and perfectly smooth walls. It is 
compared to the glass transition temperatures T s (h) of three 
studies: (i) Monte Carlo simulations of a lattice model for 
free-standing atactic polypropylene (PP) films [75| (ii) Exper- 
iments of supported atactic polystyrene (P S) fil ms (spin cast 
from toluene solution onto silicon wafers) [l28l ] (ii i) Experi - 
ments of supported, high-molecular weight PS films [I29lll30| . 
The solid line indicates T g (/i)/T g = 1 - (hq/h) s (h is a ma- 
terial dependent characteristic lengt h) |l3ll . Il32| , the dashed 
line T g (h)/T g = 1/(1 + h /h) [l2a | and the dotted line the 
approximation 1 — ho/h, valid for small ho/h. 



aspect is demonstrated in Fig. [5] for the case of the poly- 
mer model of Eqs. {1} and ([2]) close to a purely repulsive 
and perfectly smooth substrate (Eq. ((5]) with e w = e and 
fw = 0). The figure shows the structural relaxation time 
versus the distance from the substrate for various temper- 
atures. While the effect of the substrate is rather weak 
and short ranged at high temperatures, the strength of 
the substrate effects grows significantly as the tempera- 
ture is decreased. Furthermore, the spatial extension of 
the region affected by the substrate also becomes larger 
upon cooling. At the highest temperature shown, the 
substrate effects are visible only within a distance of 2-3 
particle diameters. In contrast to this, the range of sub- 
strate effects exceeds 6 particle diameters at the lowest 
temperature investigated. The most remarkable effect, 
however, is on dynamic properties: At the lowest tem- 
perature investigated, the relaxation time decreases by 
roughly two orders of magnitude when approaching the 
substrate from infinity. 

It is noteworthy that the above observations via com- 
puter simulations of a very simple model are qualita- 
tively in l ine with experimental findings on real polymers 
jl26l \l2i\ . This indicates that, at least on a qualitative 
level, generic features such as geometric confinement as 
well as adhesive / repulsive nature of the substrate play a 
more important role than material specific details. 

The above discussed effects of substrate on the dy- 
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FIG. 7: Components of the radius of gyration and the end- 
to-end distance in directions parallel (R g » and R^ e y) and 
perpendicular (R% ± and R^ ±) to the substrate versus the 
distance, z, from the substrate. R^ y (z) (left ordinate) and 
R^ e || (z) (right ordinate) behave qualitatively similar. They 
develop a maximum close to the substrate and then converge 
towards a constant (bulk) value in the film center indicated by 
horizontal dashed lines. (Note that Tig y and R gt ± are shifted 
upwards by an amount of unity in order to avoid crossing with 
the end-to-end curves) 16011. 



namics of structural relaxation have strong implications 
regarding the thermal and mechanical properties of the 
system. Indeed, in the case of perfectly smooth and re- 
pulsive walls, the confined system behaves more liquid 
like compared to the same polymer system in the bulk 
(infinitely far from the substrate) . Similarly, the presence 
of atomistically corrugated substrates with adhesive in- 
teractions increases the solid character of the system. For 
sufficiently small slab or film thickness (h < 100 nm), this 
behavior translates itself into a dependence of the glass 
transition temperature T g (the temperature at which the 
polymer forms an amorphous solid) both on the type and 
thickness of the slab. Figure[6]illustrates this issue for the 
case of smooth non adhesive walls, where the expected 
reduction in T g is observed in experiments accompanied 
by a similar reduction in T c in simulations (T c is the so 
called ideal glass transition tem perature within the mode 
coupling theory [l33l . 1 134L Il35l |). It is noteworthy that, 
in the case of adhesive substrates, the expected opposite 
effect, i.e. an increase in the glass transition temperature 
upon confinement is indeed observed (see e.g. Torres and 
coworkers [82| and references therein). 

Let us now turn to the question how polymer spe- 
cific properties such as the conformation of a chain 
may change close to a substrate. For this purpose, we 
show in Fig. [7] that polymer conformation is significantly 
stretched along the parallel direction in the proximity of 
the substrate. For this purpose, we compare the parallel 
and perpendicular components of the radius of gyration, 
R% || , JJg j_ as well as the same components of the end- 



iig e j_ in a slab geometry. 



to-end distance, i?^ e ^ , 

The components are plotted versus the distance, z, 
from the substrate, where z denotes the position of 
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chain's center of mass. So, R? Az), for instance, is the 
radius of gyration parallel to the substrate, which is av- 
eraged over all chains whose centers of mass are located 
at z. The figure shows that both the radius of gyration 
and the end-to-end distance agree with the bulk value 
if z > z w + 2i?g Ulk . Here, z w « 1 is the wall position, 
i.e. the smallest distance between a monomer and the 
substrate. As the chain's center of mass approaches the 
substrate. R~ ,, and R? n first develop a shallow mini- 
mum and then increase to about twice the bulk value 
followed by a sharp decrease to zero in the very vicinity 
of the substrate where practically no chain is present. On 
the other hand, the perpendicular components, R^ ± and 
i?g e j_ , first pass through a maximum before decreasing 
to almost at the substrate. This beha vior has been 
observed in several other simulations (see |l36l | and ref- 
erences therein) , also for larger chain length than studied 
here ^37l fl38j . 

VI. CONCLUSION 

In this paper, we provide a brief survey of modeling 
and simulation studies of dense systems of flexible lin- 
ear polymers close to solid substrates. The challenging 
character of polymer science relies upon the fact that 
the smallest and largest length scales present in a poly- 
mer system may span many orders of magnitude, usually 
~ 1A (size of an atom) up to hundreds of nanometers 
(end-to-end distance) . This broad range of length scales 
brings about a correspondingly wide time window. An 
adequate study of polymer systems thus necessarily in- 



volves the use of simplifying concepts allowing one to 
focus on essential features. It is, therefore, not surpris- 
ing that coarse graining procedures are common practice 
in polymer science. We therefore address in this arti- 
cle some of the basic ideas regarding the development of 
coarse grained models for a computational study of poly- 
mer systems. Along this route, we also address Monte 
Carlo methods as compared to molecular dynamics sim- 
ulations. 

Some salient properties of polymers close to a solid 
substrate are also presented. In particular, we show how 
the presence of a solid substrate may affect both the 
static and dynamic properties of a polymer melt. This 
is exemplified by a significant slowing down of the local 
diffusion dynamics -and the closely related dynamics of 
structural relaxation- in the proximity of attractive sub- 
strates. Similarly, a generic enhancement of diffusion (re- 
sulting in a reduction of the glass transition temperature 
in a slab geometry) is observed close to perfectly smooth 
non adhesive surfaces. These findings are evidenced both 
experimentally and by computer simulations. Polymer 
specific features, on the other hand, reflect themselves, 
e.g. in a change of conformational degrees of freedom 
from fully isotropic in the bulk to an elongated (in direc- 
tion parallel to the substrate) state in the very vicinity 
of the substrate. 

We gratefully thank J. Baschnagcl, Wolfgang Paul and 
Marcus Miillcr for useful discussions and providing us 
with results of their recent work. FV acknowledges the 
support from the industrial sponsors of the ICAMS, the 
state of North Rhine- Westphalia and European Union. 
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